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Abstract 

We calculate the Meissner masses of gluons in neutral three-flavor color superconducting matter 
for finite strange quark mass. In the CFL phase the Meissner masses are slowly varying function of 
the strange quark mass. For large strange quark mass, in the so called gCFL phase, the Meissner 
masses of gluons with colors a = 1,2,3 and 8 become imaginary, indicating an instability. 

PACS numbers: 12.38-t 
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I. INTRODUCTION 



At asymptotic densities cold quark matter is in the Color Flavor Locked (CFL) phase 
of QCD [lj (see also Q and 0). This state is characterized by nine gapped fermionic 
quasi-particles (3x3, for color and flavor) and by electric neutrality even for non vanishing 
quark masses Mj, provided Mj j£ does not destroy the CFL phase j^J. For lower densities, 
it has recently been shown P, Q, Q] that, including the strange quark mass M s , requiring 
electrical and color neutrality, and imposing weak equilibrium, a phase transition occurs, 
from the CFL phase to a new phase, called gapless CFL or gCFL. In the gCFL phase 
seven fermionic quasiparticles have a gap in the dispersion law, but the remaining two are 
gapless. At zero temperature the transition from the CFL to the gCFL phase takes place 
at Mg/fij, ~ 2 A, where yn> is the quark chemical potential and A is the gap parameter. At 
non zero temperature the situation is more involved M, and also the existence of mixed 

n 

phases 191 has to be taken into account. The next phase at still lower densities is difficult to 

n 

determine and the crystalline color superconductive phase is a candidate [1Q |. 

The aim of this paper is to investigate the dependence of Meissner masses on the strange 



quark mass in the gCFL p 
by Huang and Shovkovy 



lase. For two flavors a similar analysis has recently been performed 



111 ]. Imposing weak equilibrium and neutrality they compute the 



gluon Meissner mass in the 2SC phase and show that an instability arises in a certain 
range of values of the parameters, with some gluon masses becoming imaginary. We present 
numerical evidence that a similar instability is also present in the gCFL phase. By the same 
method we also investigate the dependence of Meissner masses on M s in the gapped phase 
CFL with M s ^ 0). Our calculational scheme is the High Density Effective Theory (HDET) 



131 ]. which allows a significant reduction of the computational complexity. 



The plan of this paper is as follows. In Section |H] we derive the effective Lagrangian 
and the grand potential in the HDET scheme. In Section ITTT1 we determine the polarization 
tensor for the gluons in the HDET approximation. In Section IIVI we present and discuss 
the numerical results for the Meissner masses as a function of the strange quark mass. The 
conclusions are summarized in Section 
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II. HDET APPROACH TO THE GCFL PHASE 



Following Ref. p the Lagrangian for gluons and ungapped quarks with M u = Ma = 
and M s ^ can be written as follows (color, flavor and spin indices suppressed): 

C = i) (z^-M + z^o) ^ (1) 

where M = diag(0, CLM S ) is the mass matrix in flavor space and the matrix of chemical 
potential is given by [5| 



/iff = (ii b 8 i:j - n Q Qij) S a(3 + 5 i:j + ^=^T^\ 



(2) 



= 1,3 flavor indices; a,/3 = 1,3 color indices). Moreover T 3 = |diag(l, — 1, 0), 
T 8 = ^gdiag(l, 1, — 2) in color space and Q = diag(2/3, — 1/3, — 1/3) in flavor space; /iq 
is the electrostatic chemical potential; /i 3 ,/i 8 are the color chemical potentials associated 
respectively to the color charges T 3 and T 8 ; ^ is quark chemical potential which we fix to 
500 MeV. As usual in the HDET, to get rid of the Dirac structure we introduce velocity 
dependent fields of positive (negative) energy -?/> v (\I/ v ) by the Fourier decomposition j3] 

^E e<WV ' X (V>v + *v) , (3) 

V 

where v is a unit vector representing the Fermi velocity of the quarks. Substituting the 
expression Q in the Eq.(Q) at the leading order in M^/fib we obtain the HDET Lagrangian 



£ = X>t \iV-D + 5fi- — U v - P^t 
where W = (1, v), = (1, -v) and 



D»D V 



JV-D + 2fi t 



A, (4) 



P v» = g ^ -- (V»V U + WV") . (5) 
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It is clear that, at this order of approximation, the effect of M s ^ is to reduce the chemical 
potential of the strange quarks. Let us now define a new basis ipA for the spinor fields: 

v«< = E ^) ai i>a , (6) 

A=l 

where the matrices F A can be expressed by 

Fx = -I + T 3 + -^=T B , F 2 = -I - T 3 + -^=T 8 , F 3 = -I -^=T 8 , (7) 
3 V3 3 ° V3 3 ^ 

F 4 , 5 = Tx±iT 2 , F 6J = T 4 ±iT 5 , F 8i9 = T 6 ±zT 7 , (8) 



with T a = A a /2 the SU(3) generators and Iq the identiy matrix. Introducing the Nambu 
Gorkov fields 

1 I V 



Xa 



V2 \ C^*_, 
the kinetic part of the Lagrangian (j3J) reads 



(9) 



£o = E Xa 



1 f V.l+6p A -$)fi. w 











2fi b 



>AB 



XB 



where 



5 ha = (Sh™, S/jgd, Sfx bs , 5/j, rd , 5/i du , 8/i rs , 5/j bu , Sfx gs , Sum) 



and M\ = M s 2 (0, 0, 1, 0, 0, 1, 0, 1, 0). If we define 



5 ha = 



Ml 
2/i 6 



(10) 



(11) 



(12) 



we may recast Eq. (fTUj) as 



£o = E Xa 



( 



(V-1+ 5hT) Sab 



V 







Xb 



(13) 



(V-£- 5nf) Sab 

which is formally equivalent to the Lagrangian for massless quarks with different chemical 
potentials. 

In the gapless color-flavor-locking (gCFL) phase the symmetry breaking is induced by 
the condensate 

A°f =< Via C 75 Vffi >= E A ' ^ ^ ( 14 ) 



1=1 



and the corresponding gap term in the Lagrangian in the mean field approximation is 

£a = -- A°f E €,-v C 75^-ft+v + h.c. . 

Z v 

In the Nambu-Gorkov basis (jHJ) the gap term reads 



(15) 



£a = E Xa 



where Aab is the 9x9 matrix defined by 



-A 



AB 



-A 



Xb 



(16) 



AB 



A AB = -E A/ Tr [fJ ej F B ej 
i=i 
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(17) 



From Eqs. ljlHj) and (fTKj) one immediately obtains the inverse fermion propagator that in 
momentum space is given by 



SabW 



(V-£ + 8fi e J f )5, 



Mi —A AB 

eff\ 



-A AB (y-l-8nT)8AB 
It can be inverted to get the fermion propagator 



(18) 



WO 



(PA(V-l-5fi eff ) A- 1 
Dac Acb 



where 



AB 



1 



Pac A C b 



{DA(V -l + Sfi 



eff 



A" 1 ' 



AB 



(19) 



(20) 



A (V ■ I - 5fJL e f f ) A" 1 (V ■ I + 5fl e ff) - A 2 

and D = P(V <-> V,5fi e ff <-> —Sfi e ff). From the poles of the propagator we can now 
determine the dispersion laws of the quasiparticles. The knowledge of the dispersion laws 
allows the evaluation of the grand potential which in the limit of zero temperature is given 
by 



n 



i 



2tt 2 Jo 



3=1 



i/7\ 2 4- A 2 + A 2 ^i - ^ Q 



(21) 



where A is the ultraviolet cutoff, ej ) are the quasi particle dispersion laws, i\\ is the quark 
momentum measured from the Fermi surface (p = fib + and G is the Nambu-Jona 
Lasinio coupling constant. We fix G as in |5j, requiring that the value of the gap is 25 Mev 
for M s = 0. In order to enforce electrical and color neutrality one has to minimize the grand 
potential with respect to fiq, fi 3 and fi 8 . Including the stationary conditions with respect to 
the gap parameters Ai, A 2 , A 3 one ends up with a system of six equations which must be 
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solved simultaneously. Once this system of equations is solved one may express Ai, A 2 , A 3 
and the chemical potentials [Iq, ^3 and fis as a function of M^/fib- We have numerically 
checked that using the grand potential (j^D), with fi b = 500 MeV and A = 800 MeV, we 
recover the results of Ref. 51 with an error of 5%. 



III. POLARIZATION TENSOR OF GLUONS 



To compute gluon Meissner masses we evaluate the polarization tensor H^{p). In the 
HDET approach, at the leading order in gfif,, there are two contributions to the polarization 
tensor: The self-energy (s.e.) diagram and the tadpole (tad) diagram (see e.g. Fig. 2 in 
To evaluate the self-energy diagram we extract the trilinear quark-gluon coupling by 
the minimal coupling term in the Lagrangian (0J): 



V 

which can be rewritten in the Nambu-Gorkov basis (jHJ) as 



(22) 



A = ig E x\ 



Xb = ig^2 Xa H'ab Xb A a (23) 



f i V ■ A a h a AB 
v -iV-A a h A * B 

where h a AB = Tr[F A T a F B ]. Therefore the self-energy contribution to the polarization tensor 
is given by: 



in 



s.e.fiu 



(P) 



g Vt f rfv 



Att a J An 



— / d z lTr 



S{t) H afl S{£ + p) H 



bv 



(24) 



In order to evaluate the tadpole diagram contribution we extract the quadrilinear quark- 
gluon coupling from the second term on the r.h. side of Eq.(jIJ) 



r 2 ; t T a Tb 

C-2 = -g Z^WvTy 



V-£ + 2/i b 

In the Nambu-Gorkov basis this term reads 



Vv Put, A^ a A v b . 



(25) 



C2 = -g 2 E x\ 



/ A Ft 



V-£+2fi b 









JaS>* 
AB 



Xb iV K K = -g 2 E X \ Yab Xb P,u A* Al , (26) 



V-i+2/j, b 



with d AB = Tr[F A T a T b F B \. The tadpole contribution is then evaluated to be 



tad,/iu 



-2g 2 



An r dv 



16tt 4 7 4tt 



P"" / dldh Tr 



( 



iS(l) 







Y 



ab 



(27) 
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Finally, in the HDET approximation, the gluon polarization tensor is given by: 

KM = nx^(p) + nif^(p). (28) 

We note that this result for the polarization tensor is correct at the order OiM s j fib) 2 ■ 

IV. NUMERICAL RESULTS 

The Meissner masses of the gluons are obtained by the eigenvalues of the polarization 
tensor (|28|) in the static limit po = 0, p — > 0. In the CFL phase with M s = the Meissner 



masses are degenerate and one has 



16] 



"4 = ^¥f----ln2 + i). (29) 
7T 2 V 36 27 2 ' 



For a non zero strange quark mass the integrals in Eqs. and (}2*Tj) have to be evaluated 
numerically In Fig. Q we present the results for the squared Meissner masses of gluons 
with color a = 1, 2, 3, 8 in units of m 2 M . The solid line denotes gluons with color a = 1, 2; 
the dashed line gluons with color a = 3; finally the dot-dashed line gluons with a = 8 (the 
physical eighth gluon is obtained by a mixing with the photon \v\ . w\ and its mass is only 
proportional to the mass of the unrotated gluon; however in the ratio m 2 M (M s ) / 'm^(O) the 
proportionality constant cancels out). We find that increasing M 2 / ^ the degeneracy in 
the Meissner masses is partially removed. Moreover there is a discontinuity of the squared 
Meissner mass of gluons of colors a = 1,2,3,8 0, 20], which at the onset of the gCFL 
phase, i.e. for M 2 / fif, ~ 2 A, drop to negative values. Thus we find an instability in the 
gCFL phase analogous to the one observed by Huang and Shovkovy [11] in the g2SC case. 

In Fig. |21 we present the results for the gluons with color a = 4, 5 (solid line) and 
color a = 6, 7 (dashed line). Also in these cases the squared Meissner mass of gluons are 
continuous functions of M 2 /fib- One can notice that, for very large values of the strange 
quark mass, the squared Meissner masses of these gluons become negative. However this 
result is not robust because in the computation of the polarization tensor Eq. (j28|) we have 
discarded terms of order 0(M S / fi b ) 2 . Therefore, to establish the instability related to the 
gluons a = 4, 5, 6 and 7 a more accurate analysis would be needed. Also in this with 
previous Fig. 1, our results give not only the Meissner masses in the gCFL phase, but also 
their dependence on the strange quark mass in the CFL phase, i.e. for M 2 /fi b < 2A. 
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FIG. 1: Squared values of the Meissner masses in units of m 2 M (see Eq. (|29[>) as a function of 
Ml I (in MeV) for gluons a = 1,2,3,8. The solid line denotes the gluons with colors a = 1,2. 
The dashed line denotes the gluons with color a = 3; finally, the dot-dashed line is for a = 8. 




20 40 60 80 100 120 M> b 



FIG. 2: Squared values of the Meissner masses in units of function of Mg/fMf,. Dashed 

line denotes the gluons with colors a = 4, 5; solid line the gluons with colors a = 6,7. 

The instability we have found means that the vacuum was not correctly identified. A 
possible origin of the instability is a non vanishing vacuum expectation value (vev) of one (or 
more) time components of the gluon operator < A° a >^ (see j^). Clearly defining a 
new field operator with vanishing vev's adds contributions that, for a = 3, 8, act as effective 
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chemical potential terms in the lagrangian: ~ g < A® > ip^Xaip. The presence of these 
new terms would alter the previous results and may lead to real Meissner masses. We have 
numerically checked that, either with < A® >^ and the other vev's < A° a >= 0, or with 
< >^ and the other vev's equal to zero, one removes the instability (numerically the 
non vanishing vev's must be of the order of ~ 10 MeV). At present the physical mechanism 
at the basis of this gluon condensation is still unclear and we do not push the analysis any 
further since our purpose here is to indicate the instability in gCFL and not to fully discuss 
its antidotes (for possibly relevant discussion see 22j). m an y case > given the instability 
of the gCFL phase, other patterns of condensation, e.g. spin-one color superconductivity, 
should be also considered (for a recent analysis see 23j and references therein). 



V. CONCLUSIONS 



It is well established that at asymptotic large densities quark matter is in the CFL phase. 
At lower densities, in a range presumably more relevant for the study of compact stars, 
neutrality, together with finite strange quark mass, suggests the gCFL (gapless CFL) phase 
as the next occurring ground state. Our calculations in this note suggest an instability of 
the gCFL phase, a phenomenon analogous to what observed in the two flavor case. The 
instability arises because gluons of color indices 1, 2, 3 and 8 present an imaginary mass. 
Its removal may require a different condensation pattern, most probably including gluon 
condensation. 
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